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Many classical scalar field theories possess remarkable solutions: coherently oscillating, localized 
clumps, known as oscillons. In many cases, the decay rate of classical small amplitude oscillons is 
known to be exponentially suppressed and so they are extremely long lived. In this work we compute 
the decay rate of quantized oscillons. We find it to be a power law in the amplitude and couplings 
of the theory. Therefore, the quantum decay rate is very different to the classical decay rate and 
is often dominant. We show that essentially all oscillons eventually decay by producing outgoing 
radiation. In single field theories the outgoing radiation has typically linear growth, while if the 
oscillon is coupled to other bosons the outgoing radiation can have exponential growth. The latter 
is a form of parametric resonance: explosive energy transfer from a localized clump into daughter 
fields. This may lead to interesting phenomenology in the early universe. Our results are obtained 
from a perturbative analysis, a non-perturbative Floquet analysis, and numerics. 



I. INTRODUCTION 

Although there is no direct evidence yet, there are good 
reasons to think that scalar fields are plentiful in nature, 
such as the Higgs boson, axion, inflaton, p/reheating 
fields, moduli, squarks, sleptons, etc. In recent years, 
there has been increasing interest in certain kinds of long 
lived structures that scalar fields can support. In particu- 
lar, under fairly broad conditions, if a massive scalar field 
<f) possesses a nonlinear self-interacting potential, such as 
4> 3 + . . . or — 4> A + . . ., then it can support coherently oscil- 
lating localized clumps, known as oscillons. Due to their 
oscillations in time and localization in space, they are 
time-dependent solitons. 

Remarkably, oscillons have a long lifetime, i.e., they 
may live for very many oscillations, despite the absence 
of any internal conserved charges. Related objects are 
Q-balls pP, which do carry a conserved charge. If the 
field <f> is promoted to a complex scalar carrying a U(l) 
global symmetry (f> — > e za (f>, then oscillons correspond to 
oscillatory radial motion in the complex <f> plane, while 
Q-balls correspond to circular motion in the complex <f>- 
plane. Since oscillons do not require the U(l) symmetry 
for their existence (in fact they are typically comprised 
of only a real scalar, as we will assume), then oscillons 
appear even more generically than Q-balls. 

As classical solutions of nonlinear equations of motion, 
oscillons (related to "quasi-breathers" ) possess an asymp- 
totic expansion which is exactly periodic in time and lo- 
calized in space, characterized by a small parameter e. 
In a seminal paper Segur and Kruskal [2] showed that 
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the asymptotic expansion is not an exact solution of the 
equations of motion and in fact misses an exponentially 
small radiating tail, which they computed. They showed 
that oscillons decay into outgoing radiation at an expo- 
nentially suppressed rate ~ exp(— b/e), where b = 0(1) is 
model dependent. Although their work was in 1 spatial 
dimension, similar results have been obtained in 2 and 
3 spatial dimensions [3]. Furthermore, various long lived 
oscillons have been found in different contexts, including 
the standard model @HS], abelian- Higgs models [7J, ax- 
ion models [5] , in an expanding universe [3 [TU] , during 
phase transitions [TTHT4] , domain walls [15] , gravitational 
systems [TMTB] (called "oscillatons"), large +(f> 6 models 
[19] . in other dimensions [^Hj, etc, showing that oscillons 
are generic and robust. 

In almost all cases, investigations into oscillons have 
thus far been at the classical level. There is a good rea- 
son for this approximation: the mass of an oscillon M osc 
is typically much greater than the mass of the individual 
particles m^. This means that the spectrum is almost 
continuous and classical. But does this imply that ev- 
ery property of the oscillon, including its lifetime, is ad- 
equately described by the classical theory? To focus the 
discussion, suppose a scalar field <f> in d + 1 dimensional 
Minkowski space-time has a potential of the form 

V&) = \m*<f + ^m^/V + ^™ 3 - d 4 + • ■ • (1) 

Since we will keep explicit track of H in this paper, it 
has to be understood that m here has units T -1 ; the 
quanta have mass — hm, and {A§, A4} have units 
hr l . It is known that for A = |A§ — A4 > such a theory 
possesses classical oscillons that are characterized by a 
small dimensionless parameter e and have a mass given 
by M OBC ~ fl(e) = /(e), where 1(e) = e 2 ~ d in the 
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standard expansion that we will describe in Section [TT| 
but can scale differently in other models (e.g., see [loj ). 
So if A h < 1(e), then M oac > m . 

In such a regime we expect that various properties of 
the oscillon are well described classically, such as its size 
and shape. In this work we examine whether the same is 
true for the oscillon lifetime. We find that although the 
oscillon lifetime is exponentially long lived classically, it 
has a power law lifetime in the quantum theory controlled 
by the "effective H" - in the above case this is A3 h or A4 h. 

In this paper we treat the oscillon as a classical space- 
time dependent background, as defined by the e expan- 
sion, and quantize field / s in this background. We go to 
leading order in h in the quantum theory. We find that 
oscillon's decay through the emission of radiation with 
wavenumbers k — 0(m). We show that the decay is a 
power law in e and the couplings, and we explain why 
this is exponentially suppressed classically. Our analysis 
is done for both single field theories, where the emitted 
radiation typically grows at a linear rate corresponding 
to 3 — > 2 or 4 </> — > 2cf) annihilation processes, and 
for multi-field theories, where the emitted radiation often 
grows at an exponential rate corresponding to <j) — > 2 \ 
decay or 2 cf) — >• 2 \ annihilation processess. We calculate 
the quantum decay rates in several models, which is sup- 
ported by numerical investigations, but our work is also 
qualitative and of general validity. We also comment on 
collapse instabilities for k — O(em) modes, whose exis- 
tence is model dependent. 

The outline of this paper is as follows: In Section |TT] 
we start with a review of classical oscillons and describe 
their exponentially suppressed decay in Section |III| In 
Section |IV| we outline the semi-classical quantization of 
oscillons and derive the decay rate of oscillons in Section 
|y] Having started with single field models, we move on to 
examine the effects of coupling to other fields in Section 
|VI| In Section |VII| we discuss when the decay products 
grow linearly in time and when it is exponential. Here we 
demonstrate that coupled fields can achieve (depending 
on parameters) explosive energy transfer, which may have 
some cosmological relevance. We comment on collapse 
instabilities in Section IVIIII and conclude in Section llXl 
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II. CLASSICAL OSCILLONS 

In this section we review the asymptotic expansion for 
a single field oscillon, valid in 1, 2, and 3 dimensions. 
We then explain why the oscillon slowly radiates at an 
exponentially suppressed rate. 

Consider a single scalar field (f> in d + 1 dimensional 
Minkowski space-time (signature H . . .) with classical 



S = d d+ 



2 ±2 

-m <p 



Vi{4>) 
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where Vi(<j)) is a nonlinear interaction potential. Here m 
has units T~ 1 . For simplicity we measure time in units of 
1/m, so without loss of generality, we set m = 1 from now 
on in the paper, unless otherwise stated. The classical 
equation of motion is 



0- V 2 + + V/(0) = 0. 



(3) 



In this paper we will focus on the following types of in- 
teraction potentials 



A3 
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where A = |Ag — A4 > is assumed. This includes the 
cases (i) A3 > and A4 > 0, which will occur for a generic 
symmetry breaking potential, and (ii) A3 = and A4 < 0, 
which is relevant to examples such as the axion. In case 
(ii) higher order terms in the potential, such as +<fi 6 are 
needed for stabilization. Note that in both cases the in- 
teraction term causes the total potential to be reduced 
from the pure </> 2 parabola. This occurs for <ft < in the 
(j) 3 case, and for |0| > in the (p A case. It is straightfor- 
ward to check that a ball rolling in such a potential will 
oscillate at a frequency lower than m = 1. If the field <j) 
has spatial structure, then one can imagine a situation in 
which the gradient term in eq. ^ balances the nonlinear 
terms, and a localized structure oscillates at such a low 
frequency - this is the oscillon. 

Naively, this low frequency oscillation cannot couple 
to normal dispersive modes in the system and should be 
stable. Higher harmonics are generated by the nonlin- 
earity, but resonances can be cancelled order by order in 
a small amplitude expansion. Let us briefly explain how 
this works - more detailed descriptions can be found in 
Refs. [HJH2]. 

In order to find periodic and spatially localized solu- 
tions, it is useful to rescale time t and lengths x = |x| 
to 



t VT 



P = xe, 



(5) 



with 0<t<la small dimensionless parameter. Here we 
search for spherically symmetric solutions. The equation 
of motion |3]) becomes 

(1 - e 2 )9 TT - e 2 (d pp <f> + d ~^d p ^j + + V/(0) = (6) 

To obtain an oscillon solution, the field <j) is expanded as 
an asymptotic series in powers of e as 



bUP,r) = J2e n cj )n (p,T). 



(7) 
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FIG. 1: Leading order oscillon profile /(p) in d = 1 (blue), 
d — 2 (red), and d = 3 (green) dimensions. We allow p to 
take on both positive and negative values here, i.e., a 1-d slice 
through the origin. 

The set of n depends on Vi{4>). For V} ~ A 3 3 ; n — 

1, 2, 3, . . ., and for Vj ~ — | A4 1 4 ; n — 1, 3, 5, Upon 

substitution of the series into eq. Q, the leading term 
must satisfy 

drr^l +<f>l=0, (8) 

with solution (pi — /(p)cosr, where f{p) is some spatial 
profile. Since r = t \Jl — e 2 the fundamental frequency 
of oscillation is evidently w = \fl — e 2 < 1. 

The next order terms in the expansion must not be 
resonantly driven by <f>\ , or the solution would not be pe- 
riodic. By writing down the equations for (f>2 and </>3 and 
demanding that the driving terms are non- resonant, we 
establish an ODE for f(p). Extracting the A dependence 
by defining f(p) = Af{p)/\/\, the ODE is found to be 

d PP f + ^ d pf - f + 2 / 3 = °- ( 9 ) 

This ODE possesses a localized solution for d = 1,2,3. 
In d = 1 the solution is known analytically f(p) — sech p, 
but is only known numerically in d — 2 and d = 3. (For 
d = 2, 3 there are infinitely many solutions, of which we 
take the fundamental solution). We plot f(p) in Figure 
[T] Altogether this gives the leading order term for the 
oscillon 

oac (p,r) = -JLf( p )coST + J2 £n <t>n(p,r), (10) 

V TI>1 

where higher order terms can be obtained in a similar 
manner. At any finite order in e this provides a peri- 
odic and spatially localized oscillon. By integrating the 



oscillon's energy density over all space, the total mass 
of the oscillon is found to scale to leading order in e as 

M oac ~ me 2 - d /\. 



III. CLASSICAL RADIATION 

In 1987 Segur and Kruskal [2 found that the above 
asymptotic expansion is not exact, as it misses an expo- 
nentially small radiating tail. They computed this using 
matched expansions between the inner core of the oscil- 
lon and infinity, involving some detailed analysis. Here 
we describe the physical origin of this classical radiation 
in simple terms. 

The existence of outgoing radiation is ultimately tied 
to the fact that the oscillon expansion in not an exact 
solution of the equations of motion, it is only an asymp- 
totic expansion which is correct order by order in e, but 
not beyond all orders. Consider the oscillon expansion, 
truncated to order N, i.e., 

N 

</> osc (x,i) = 5>"0 n (x,i). (11) 

n=l 

Lets substitute this back into the equation of motion. 
This will involve many harmonics. By construction the 
fundamental mode cost = cos(ujt) will cancel, since the 
spatial structure was organized in order to avoid such a 
resonance. However, we do obtain remainders including 
higher harmonics. This leads to the following remainder 

J(x,t) = j(x) cos(nuj t) + ... , (12) 

where for d — 1 we have 

j(x) = C N e N+2 sech N+2 (xe) + ..., (13) 

where we have included only the coefficient of the next 
harmonic, with some coefficient Cjv- n — 2 for asymmet- 
ric potentials and n — 3 for symmetric potentials. 

Let us now decompose a solution of the equations of 
motion </> aol into an oscillon piece and a correction S 

c/) BOl (x,t) = <j> 0!ic (x,t)+S(x,t) (14) 

with S taking us from the oscillon expansion to the nearest 
solution. We find that the remainder J(x,t) acts as a 
source for the correction S 

6 -V 2 5 + 5 = -J{x,t) (15) 

where we have ignored non-linear terms (5 2 etc) and a 
parametric driving term V" '(4> osc )5 (although such a term 
will be important in the quantum theory). The solution 
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in the far distance regime is obtained by standard means 
d d kduj J(fc,w)e l ( k - x -<^ 



8{x,t) 



(2 7 r) d + 1 -cu 2 + k 2 + 1 ± ^.0^ 
cos(fc rad ir + 7) 



(16) 



x (d-l)/2 



■ cos(cj rad i)j(fc rad ), (17) 



as x — > 00 (7 is a phase) , suggesting that there is a radi- 
ation tail with an amplitude determined by the Fourier 
transform of the spatial structure of the source, evaluated 
at fc rad . The kinematics of ( 16 ) says that the radiation tail 



has frequency and wavenumber: 

iu ra d =nui, fc rad = y n 2 uj 2 - 1, 



(18) 



where oj = yl — e 2 is the fundamental frequency of the 
oscillon. This is the so-called "quasi-breather" which is 
used to match onto a radiating oscillon [5] . 
Let us now evaluate j'(fc). In d = 1 wc find 



C N k N+1 sech 



7T k 

~2~e 



(19) 



So for k — 0(e) the remainder is a power law in e, but 
for k = 0(1), the prefactor is 0(1) with a sech function 
evaluated in its tail. For k — 0(1), which is the case for 
k = fc rad ; we note that the source is comparable to the 
Fourier transform of the oscillon itself, i.e., 

m - <m*o (20) 

Lets now compute the spatial Fourier transform of the 
oscillon. Taking only the leading order piece from eq. ( 10 ) 
at t = we have 



= (fc) 



4e 

71 



d d x f(xe)e 



ik-x 



(21) 



For d = 1 we have f{p) = sech p, allowing us to compute 
the Fourier transform analytically: 



= (fc) 



47T 

71 



sech 



IT k 
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(22) 



For fc 3> e this is exponentially small. The same is true 
in d = 2 and c? = 3 which can be computed numerically. 
We can summarize the behavior for d = 1,2,3 for k 3> e 
by the following scaling 



»(*) 



1 



^(efc)^- 1 )/ 1 



■ exp 



Crffc 

e 



(23) 



where the constant in the argument of the exponential 
is dimensional dependent: c\ = 7f/2, C2 ~ 1.1, and C3 ~ 
O.6. 1 If we evaluate this at the relevant mode fc = fc rad s» 



1 The constant is in fact the first simple pole of /(p) along the 
imaginary p-axis; see Ref. [3] for comparison. 



\Jn 2 — 1, the Fourier amplitude is exponentially small as 
e^O. 

Hence the longevity of the oscillon is due to the fact 
that the spatial structure is dominated by small fc = 
0(e) wavenumbers, while any outgoing radiation occurs 
through fc = 0(1) wavenumbers, whose amplitude is sup- 
pressed. The rate at which energy is lost scales as 



dE n , 



dt 



,c(fc iad )| 5 



Ac 



1 

d3T ex P 



(24) 



where b = 2 Cd Vn 2 — 1 is an 0(1) number. 

Now the energy of an oscillon scales as E oec ~ 
l/(\e d ~ 2 ), so the decay rate T d — \E~j dE osc / 'dt\ scales 
as Yd ~ e _1 exp(— b/e). The constant of proportionality 
in is non-trivial to obtain. The reason is the following: 
at any finite order in the e expansion, the corresponding 
terms in the equations of motion are made to vanish due 
to the spatial structure of the solution. However, there is 
always a residual piece left over which is 0(e°) in fc-space, 
multiplied by (p osc (k). The limit of this residual piece as 
we go to higher and higher order encodes the constant, 
but we do not pursue that here. Alternate methods for 
obtaining this constant can be found in the literature, 
such as 00 [23H2S], (see also [H 123). A special case 
is Sine-Gordon in 1-d where the remainder J vanishes as 
N — > 00, so the constant vanishes in this case. 



What we have obtained in (24) is the leading e scal- 
ing: it is from n — 2 for asymmetric potentials and from 
n = 3 for symmetric potentials. The form of the radia- 



tion in eq. ( 24 1 agrees with the scaling found by Segur 
and Kruskalfiir in d — 1 and generalized to other d. Our 
methodology of computing a spatial Fourier transform 
and evaluating it a wavenumber determined by kinemat- 
ics is general and should apply to various other oscillon 
models, such as multi-field. 



IV. QUANTIZATION 

The preceding discussion explains why a classical os- 
cillon can live for an exceptionally long time. This does, 
however, require the oscillon to be placed in the right 
initial conditions for this to occur. Some investigation 
has gone into the stability/instability of classical oscil- 
lons under arbitrary initial conditions, e.g., |28|,l29j. Here 
our focus is on a sharply defined question, free from the 
ambiguity of initial conditions: If a single oscillon, as de- 
fined by the e expansion, is present - what is its lifetime ? 
There is a sharp answer in the classical theory - expo- 
nential, and now we address the question in the quantum 
theory. 

A semi-classical (leading h) description involves treat- 
ing the oscillon 4> osc = <j) osc (t,r) as a classical background 
and quantizing fields in this background. In this section, 
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only the field <fi itself is present to be quantized, but in 



Section VI we will introduce a second field x to quantize. 
Let's write 



(x,i) =<k„(a;,t) + 0(x,t), 



(25) 



where <j> is a quantum field satisfying canonical commuta- 
tion relations. At any finite order in the e expansion, </> osc 
is an exact periodic solution of the equations of motion, 
as discussed in Section [ll| Perturbing around a solution 
allows us to write down the following equation of motion 
for <f> in the Heisenberg picture 



= 0, 



(26) 



where <I>(0 O3C ) = V"(fa so ). We have neglected higher order 
terms in <j>, since we are only interested in a leading order 
h analysis; so our results only apply in the weak cou- 
pling regime, i.e., Xh <C 1, which is a necessary condition 
for the oscillon to possess a roughly classical description. 
It seems reasonable, and is the common assumption in 
the literature on Q-balls as well as p/reheating, that for 
AS < 1 higher order effects will not modify the cen- 
tral conclusions. Eq. (26 1 is the theory of a free quan- 



tum scalar field with a space-time dependent mass. The 
ground state of this theory is given by an infinite sum of 
1-loop vacuum diagrams, coming from N$ = 0,1,2,... 
insertions of the external field <&. The AT$ = diagram 
corresponds to the ordinary zero point energy of a free 
field, while the iV$ > 1 diagrams correspond to pro- 
duction of <j) quanta from the background source. We 
will see, however, that although these are 1-loop effects 
in the semi-classical Lorentz violating theory, these pro- 
cesses will have an interpretation in terms of tree-level 
effects of the underlying Lorentz invariant theory. 

For small e, the oscillon is wide, with fc-modes concen- 
trated around k = 0(e) <C 1, which suggests it is more 
convenient to perform the analysis in /c-space. So let's 
take the Fourier transform: 



fa + U> k fa + 



d d k' 
W) 



^(k-k')^' =0, 



(27) 



where ui k = k 2 + l. Here we used the convolution theorem 
on the final term. Now if the background was homoge- 
neous, then each fa would decouple. This makes the 
solution rather straightforward, as is the situation during 
cosmological inflation, for instance. In that case each fa 
is proportional to a single time independent annihilation 
operator d k times a mode function v k (t) that satisfies the 
classical equation of motion, plus hermitian conjugate. 
But due to the inhomogeneity in </> osc , the k- modes are 
coupled, and the solution in the background of an oscil- 
lon is non-trivial. 



Nevertheless a formal solution of the Heisenberg equa- 
tions of motion can be obtained. The key is to integrate 
over all annihilation operators: 



fa(t) = Vh 



d d q 

j^a q v qk (t) + h.c. 



(28) 



Upon substitution, each mode function v qk (t) must sat- 
isfy the classical equation of motion 



d d k' 



v qk +uj k v qk + J j^Y^ik - k')v qk , = 0. 



(29) 



We see that we have a matrix of time dependent mode 
functions v qk (t) to solve for. We choose initial conditions 
such that 4> is initially in its unperturbed vacuum state. 
This requires the following initial values of the mode func- 
tions: 



Mo) 



i 



(27r) d £ d (q-k), 



V2wfe 

v q k(0) = -iu k v qk (0) 



(30) 
(31) 



The local energy density u and total energy E in <p at 
time t can be defined by the unperturbed Hamiltonian, 



giving 



h f _d d ^d d ^d d l^_ i(k '_ k) . x 
«UMJ - 2 J (27r)d (27r)d e 



E{t) 



X [VqkV* qk > +ul k <V qk V* k 

h f d d q d d k n . |2 2 1 2 



{2ir) d (2tt) 



(32) 

-[\v qk f +UJ l\v qk \ 2 ] (33) 



(uJ kk , = k • k' + 1), where the second equation is easily 
obtained from the first by integrating over x. Initially 
the total energy is E(0) — f d d k\hbJ k S d (0); the usual 
(infinite) zero point energy of the vacuum. The energy 
corresponding to <p production is contained in the time 
evolution of E(t). As we will explain, radiation comes 
from specific wavenumbers that are 0(1), and are con- 
nected to tree-level processes of the underlying micro- 
physical theory. This makes it straightforward to extract 
the correct finite result for the produced radiation energy, 
despite the zero point energy being UV divergent. 

To solve the system numerically we operate in a box of 
volume V — L d and discretize the system as follows 



d d k' 
(27r) d 



5 d (q-k) -> 



-y 



k=^,n, : eZj. (3-1) 



V 



(2vr 



(35) 



The discretized equations of motion represent an infinite 
set of coupled oscillators. They have a periodic mass 
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as driven by the background oscillon and as such are 
amenable to a generalized Floquet analysis of coupled 
oscillators. We use the Floquet theory to solve for the 
late time behavior, as discussed in Section [VII A| 

It is interesting to compare this to a classical stability 
analysis. Here one should explore all initial conditions 
that span the complete space of perturbations. To do 
this, we can write S<pk — > v q h where q is merely an in- 
dex that specifies the choice of initial condition. To span 
all initial conditions, choose v q k(0) oc <5(q — k), as in the 
standard Floquet theory. But this is precisely what we 
have done to solve the quantum problem. Hence a clas- 
sical stability analysis over all initial conditions involves 
the same computation as solving the quantum problem 
for a fixed initial condition - the ground state. 



V. QUANTUM RADIATION 

Perhaps the most interesting oscillons are those that 
are stable against instabilities that would appear in a 
classical simulation; such instabilities typically arise from 
k = 0(e) modes of the oscillon, and will be discussed in 
Section |VIII| Classically, these oscillons appear to be ex- 
tremely stable. In this Section we calculate and explain 
why the oscillon lifetime is shortened in the quantum the- 
ory due to k = 0(1) modes, depending on the size of the 
effective h. 

In order to make progress, we will solve the mode func- 
tions perturbatively. Although our methodology is gen- 
eral, we will demonstrate this with a model that makes 
the computation the easiest. Lets consider the potential 
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(36) 



The -A</> 4 term ensures that classical oscillons exist in 



the form given earlier (eq. ( 10 ) ) . The A5 (f>° will be quite 



important in the quantum decay, but it does not effect 
the leading order term in the classical oscillon expansion. 
The dots represent higher order terms, such as +<fi 6 which 
ensure that the potential is well behaved at large (p. It is 
important to note that this potential, as well as all models 
studied in this paper, possess oscillon solutions which are 
highly stable when studied classically. The only regime 
in which the classical oscillons are unstable depends sen- 
sitively on dimensionality and amplitude, which we will 
mention in Section rVlID 

It can be verified that for this potential the oscillon has 
the form 



osc (a;, t) = (j) e (x) cos(wt) + (D(t 



where 4><l{x) = 4ef(xe)/^/\ (note 4> e ~ e). A5 only enters 
the expansion at 0(e 4 ) (generating even harmonics). 



Lets perform an expansion in powers of A5, i.e., 



5 "q 



(38) 



(i) 

with each term v q implicitly an expansion in powers of 
e, with e assumed small. We have suppressed the second 
index on v q , which would be k in momentum space or 
x in position space, since both representations will be 
informative. In position space at zeroth order in A5, we 
have 



;(0) 



72„(o) , „(o) 



A 



cj>l{x) cos 2 {wt) + 0( 



,(0) 



(39) 



The first term on the RHS is sometimes responsible for 



collapse instabilities, as we shall discuss in Section VIII 
but only for wavenumbers k = 0(e). The existence of 
such instabilities are highly model dependent and are not 
the focus of this section. Instead here we focus on k — 
0(1) which is relevant for outgoing radiation. The second 
term on the RHS, which is 0(e 4 ), is indeed relevant to 
the production of radiation, but it will be superseded 
by radiation at 0(e 3 ) that will enter when we exa mine 
v q x \ For now, we only need to conclude from eq. (39 1 
that v q ^ is equal to the solution of the free-theory plus 
0(e 2 ) corrections. In fc-space this means we take the 
unperturbed mode functions 



v 



(0; 



(*) 



(2^) d <5 d (q-k) + 0(e 2 ), (40) 



qk vv V^k~ 

which match the initial conditions mentioned earlier in 



eqs. (30 311. Note that an integral over q recovers the 



standard mode functions that occur in a free theory. 
At next order in A5 we have 



J qx 



v 2 ^ 



J qx 



^(x)cos 3 M)«(°) (41) 
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Fourier transforming to A-space and inserting v^°J gives 



V qk + UJ k V qk - 



&(q-k) 



q k 



(42) 



where ipe(k) is the Fourier transform of (f>^(x) and Q = 
3w. Here we have ignored terms on the RHS with 
frequency u) ± ui q which cannot generate a resonance. 
The first term included does generate a resonance at 
Lj q 1=3 Uk ~ uj/2. This is the equation of a forced oscillator 
whose solution can be readily obtained. By imposing the 



(37) initial conditions v ]} (0) 



J qk 



(0) = we find 



$(*) 



3!2V2^/c 



(43) 
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where 



-i(uj-\-uj q )t 



i{Gi—u q )t 



+ 



1 — — I e 



— itJkt 



-LU 2 + (U)k - W,) 2 -W 2 + (Uj k + UJ q ) 



Lets now insert our solution for v q k into the expres- 
sion for the energy density u(x, t) (eq. (32)). This in- 
cludes terms scaling as v^vty, which is the usual zero 

point energy. Another term scales as A5 v q ° k v qk , which 
is non-resonant. Then there are two important terms: 
A 2 vQvQ and A 2 



(0) (2) 
V qk V qk>- 



It can be shown that the sec- 



5 u qk u qk' 

ond term here provides the same contribution as the first 



term. (In the case of a homogeneous pump field the mode 
functions are often written in terms of Bogoliubov coef- 
ficients, which make this fact manifest. A similar argu- 
ment goes through in our more complicated inhomoge- 

neous case.) Given this, we will not write out the explicit 

(2) 

(44) solution for v qk here. Finally, we shall use the fact that 
only terms near resonance contribute significantly to the 
integral, and further we shall make simplifications using 
Lu q rs ujk ~ w/2 whenever we can; however, this must not 
be done in singular denominators or in the arguments of 
oscillating functions. 

Altogether at leading order in A5 and e we find the 
following expression for the change in the energy density 
from the zero point: 



5u(x, t) 



d d q d d k d d k' 



- d cos((k'-k).x)^(q-k)C(q-k') 



(3!) 2 2 8 J (27r) d {2ir) d (2tt) c 

1 + cos(t(cjfe/ — cuk)) — cos(i(wfc + ujq — uj)) — cos(£(wfc< + uj q — oSj) 



UJ q (0Jk 



- UJ n 



':){Uk' +0J q -U>) 



(45) 



For the 1-dimensional case we have evaluated this nu- 
merically; the results are plotted in Fig. ^ . As the figure 
shows the energy density initially grows in the core of the 
oscillon before reaching a maximum (around Su ~ 10~ 7 
in the figure for the parameters chosen). This occurs at 
each point in space and so the net effect is of continual 
growth in Su by spreading out. This has a simple inter- 
pretation: </)-particles are being produced from the core 
of the oscillon and moving outwards. We will show using 
scattering theory that this is connected to the annihila- 
tion process 3 (j> — > 2<j). 

Let us turn now to the total energy output (after sub- 
tracting the zero point) 6E(t). This comes from integrat- 
ing eq. (45) over x, giving 



SE(t) 



Agft 
(3!) 2 2 7 



d d q d d k 
(2^(2^'™" 
I - cos(i(w fe - 



•«)) 



LOq(u} k 



(46) 



This may be further simplified by recognizing that \ip t {q— 
k) 1 2 is non-negative and sharply spiked around q — k « 
for small e (as we discussed in Section III) and so it acts 
like a ^-function, i.e., 



hMq-k)| : 



d d k 
(27r) d 



l^e(k)| 2 



(27r)^ rf (q-k), (47) 



that the integration of |^ c (q — k)| 2 over all k gives the 
correct value when we replace it by the <5-function. This 
quantity has a nice interpretation when re-written as an 
integral over position space. Since ipe(k) is the Fourier 



i.x 10 



8.x 10 




where we have intoduced the integral prefactor to ensure 



FIG. 2: Energy density 5u(x,t) (in units of h). Each curve is 
at a different time interval: blue is t = 10, red is t = 25, green 
is t — 50, orange is t — 75, and cyan is t = 100. Here d = 1, 
e = 0.05, and A§/A 3 = 1. 
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transform of 4> 3 (x), we can write 



/ 



d d k 



- d \Mk)? 



d d xu 3 (x), 



(48) 



where u e (x) (f>^(x)/2 is the oscillon's local energy den- 
sity. Inserting this into ( 46 ) allows us to write the energy 



output in terms of an integral over the occupancy number 
n k (t) as follows 



with 



n k {t) 



/d d k 
j^-^huj k n k (t) 



Xljd d xu 3 (x) 1 -cos(t(2cj fc -w)) 



(3!) 2 2 4 



Lu 2 (2w k 



(49) 



(50) 



Note that the occupancy number is resonant at uj k ~ cD/2. 
The late time behavior can be studied by noting that as 



1 — cos(t(2uj k — £j)) irt 



(2ui k - a)) 2 



S(u k -u/2). (51) 



This shows an important connection; our results are a 
form of Fermi 's golden rule, where the delta- function pre- 
scription is always used and is known to provide accurate 
transition rates. 

Note that no other form of "renormalization" is re- 
quired here, since the process we are computing is con- 
nected to tree-level processes of the underlying Lorentz 
invariant quantum held theory (albeit 1-loop in the semi- 
classical Lorentz violating theory) and not loop-level pro- 
cesses which are higher order in the couplings and H. We 
will elaborate on this connection shortly. 

To compute the energy output in the Fermi golden rule 
approximation, we begin by w ritin g d d k = dfldkk d ~ 1 . 
The angular integration in ( 49 1 is trivial J dil — 
2ii d / 2 /T(d/2). Now write dk = duj k u k /k. The integral 
over duj k is trivial due to the delta-function j5l]), which 



enforces k to take on its resonant value corresponding to 
radiation A: rad , which satisfies uj k — yj fc 2 ad + m 2 = 3w/2. 
This gives the following expression for the energy output 
in the Fermi golden rule approximation 

SE(t) = j^- A . / d d x u 3 (x) t. (52) 

v > ( 3 i)2 2 4 r(|)(27r) d J eV ' y ' 

Hence the energy output increases linearly with time. 
So the oscillon must lose energy at this rate. The decay 
rate is T d = \E~l dE oa Jdt\ and using E osc — J d d xu e (x) 
we obtain 



T d (30^ 2, 



X 2 H K d / 2+1 k d ~ 2 $ d d xul{x) 
(3!) 2 2 4 r(|)(27r) d J d d xu e (x) 



(53) 



as our hnal result for the decay rate. We have labelled 
this "3 (j> —> 2 </>" annihilation for a reason we now explain. 

Its useful to connect this result to ordinary perturba- 
tion theory for a gas of incoherent particles. The differen- 
tial transition rate for a single annihilation of Ni 6 —> 2 cj> 
in a box of volume V box , evaluated on threshold for non- 
relativistic initial particles of mass m, is given by (e.g., 
see [3D]) 



yl— JVj 

<tri = „,';° x ^ \M\ 2 (2v) D s D (Pi - Pf)U 



Ni\{2m) N i 



f 



d d Pf 1 
(27r) rf 2~¥ 



/ 



■(54) 



(D = d + 1). To connect to the previous calculation 
we choose N = 3 since the relevant interaction is pro- 
vided by 3 </) —> 2 tj> annihilation at tree-level due to the 
interaction term AV7 = ^A 5 5 . This gives \M\ 2 = A 2 . 
Performing the integration over phase space (including 
division by 2 since the final 2 particles are indistinguish- 
able) gives the result 



rx(30^ 2, 



X 2 H ir d / 2 



+ljud-2 



v- 



3!2 4 r(|)(27r) d 2E 



f 



(55) 



The inverse of this is the time taken for a given triplet 
of (j)'s to annihilate. For a box of particles the total 
rate for an annihilation to occur is Ti multiplied by the 
number of indistinguishable ways we can choose a triplet, 
i.e., 



r t „ t (3(/> 



26) = 



ri(3(/>-> 2 0) (56) 



3! 



Ti(30^ 2< 



(57) 



for ^> 1 (the semi-classical regime). Since each pro- 
cess produces a pair of particles of total energy 2Et, the 
total energy output is 



6E(t) = 2E f r tot (3<f>^2 4>)t. 



(58) 



Together with (551 this recovers the result in eq. &52 



precisely, as long as we make the replacement 



v-jNi 



d d xu 3 (x). 



(59) 



Hence to leading order in A5 and e we find that decay 
rates computed using tree-level scattering theory recover 
the result computed from solving the mode functions for a 
coherent background oscillon. The important connection 
is provided by eq. ( 59 1 to account for the spatial structure 
of the oscillon. 

Given this connection, the generalization of our re- 
sult to arbitrary interactions is relatively straightforward. 
Consider the more general interaction potential 



TO) 



a 3 

3! 



A 



4! v 



A5 
5! 



Ag 
6! 



(60) 
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FIG. 3: Feynman diagrams for the process 4 0—)- 2(f). Eval- 
uated on threshold, the first diagram is iAt/4 (x 6 crossing 
symmetries), the second diagram — iAf/8 (x 4 crossing sym- 
metries), and the third diagram is — 



As mentioned earlier, a requirement for the existence of 
small amplitude oscillons is |A| — A4 > 0; this in fact is 
the only requirement on the couplings [25 . If cither A3 or 
A5 are non-zero, then the leading order radiation arises 
from 3 <p — > 2 <p annihilation, as we computed previously 
for the A3 = case. This again gives the result of eq. (53 1, 
with the generalization of the X 2 prefactor being replaced 
by the square of the matrix element for all such tree-level 
scattering diagrams, i.e., 



A? 



\M{2,cj)^ 2 0) | s 



(61) 



This includes the important case in d = 3 with A3 > 0, 
A4 > and A5 = Ag = . . . = 0, which is a renormal- 
izable field theory. For brevity, we have not drawn the 
many relevant diagrams here, but in Fig. [3] we draw the 
diagrams that are relevant to the next example. 

Now consider a purely symmetric potential with A3 = 
A 5 = and A4 < and A6 > 0. In this case, the leading 
order annihilation process is Acf) — ¥ 2 0, so we must choose 
Ni = 4 in eq. ( 54 ) . By summing the diagrams of Fig. [3] 



we find |7W(40 20)| 2 = (A| - A 6 ) 2 . Altogether we 
obtain the following result 



(A 2 



X 6 ) 2 fnr d / 2+1 ki- 2 Jd d xu 4 e (x) 



(4!) 2 2 5 r(f)(27r) d J d d xu e {x) 



(62) 

In this case the kincmatical requirement on fc rad for reso- 
nance is (jJk = \J fc 2 ad + to 2 = 2uj. Note this rate vanishes 
if and only if X§ = A|, which is true for the Sine-Gordon 
potential V SG (<p) = (1 — cos(V~X(j)))/X. In fact it can be 
proven that in (l+l)-dimensions every ip number chang- 
ing process vanishes on threshold for the Sine-Gordon 
potential. 

Lets discuss the scaling of our results ( p53| and (62 1. 
Recall that </> oac ~ e/VX, so u e ~ e 2 /A. Thus (depending 
on the relative size of couplings) the decay rates scale as 



r d (3, 
r d (4< 



->2< 
-> 2< 



e A \\n/\ 2 , 

e 6 A 4 /i, 



or e 4 Xfti/X 2 (63) 
or e 6 X 2 6 h/X 3 . (64) 



(We also find that in the large +A6 4> 6 model of Ref. [T§] 
the scaling is Td ~ X\h/ X§.) Hence the quantum mechan- 
ical decay rate is a power law in the parameters e, A;. 
Such quantum decay rates will only be smaller than the 



corresponding classical decay rate ( 24 ) if the "effective h" 



(such as A h) is extremely small, as we are comparing it 
to an exponentially small quantity in exp(— b/e). Unlike 
collapse instabilities that we will mention in Section [VIII[ 
it is almost impossible to avoid this radiation by changing 
dimensionality, parameter space, or field theory (the only 
known exception is the Sine-Gordon model in d = 1). 

One may wonder why the classical analysis failed to 
describe this decay rate accurately given that the oscil- 
lon can be full of many quanta, say « M aBC /m,p, with 
an almost continuous spectrum. Let us explain the res- 
olution. For N$ 3> 1 the quantum corrections to the 
oscillon's bulk properties are small. For example, the 
quantum correction to the oscillon's width, amplitude, 
total mass, etc, should be small, since the classical values 
are large. In fc-space we can say that these properties 
are governed by k = 0(e) modes, which carry a large 
amplitude. However, the radiation is quite different; it 
is governed by k — 0(1) modes, which are exponentially 
small in the classical oscillon, but are not exponentially 
small in the quantum oscillon due to zero point fluctu- 
ations. Another way to phrase this is to consider the 
commutation relation 

^(k)Tr(k') = 7r(k')0(k) + ih (27r) d <5 d (k - k'). (65) 

For k = 0(e) the classical value of the LHS and its coun- 
terpart on the RHS are both large, so the h correction 
is negligible. But for k = 0(1) the classical values are 
small, so the quantum corrections are important. 

An interesting issue is the behavior of the growth at 
late times. In the case of a homogeneous background 
pump field (for example, as is relevant during p/reheating 
at the end of inflation, e.g., see [3TM33| ) it is known that 
the linear growth is really the initial phase of exponen- 
tial growth. So is the same true for spatially localized 
oscillons? Here we claim that for sufficiently small am- 
plitude oscillons, the answer is no, the linear growth rate 
is correct even at late times, although this can change for 
sufficiently large amplitude oscillons where our perturba- 
tive analysis breaks down. We have confirmed the linear 
growth for small amplitude oscillons in 2 ways: (i) by 
expanding v q k to higher order in the coupling and (ii) by 
solving the full mode function equations ( |29[ ) numerically. 
The reason for this result is subtle and will be discussed 
in detail in Section [VII| we will demonstrate that whether 
the growth is in the linear regime or exponential regime 
depends critically on the oscillon's amplitude, width, and 
couplings. 
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VI. COUPLING TO OTHER FIELDS 

Most fields in nature interact considerably with others. 
It is important to know what is the fate of a </>-oscillon 
that is coupled to other fields. Let's couple (j> to another 
scalar \ and consider the following Lagrangian 



1 



„»2 i2 

77i <pr 



1 



W) + o(0x) 



1 



2 2 



1 



-gimcj)x z 



±2 2 
" 92 4> X ■ 



(66) 



The last two terms represent interactions between the 
2 fields, with gi coupling parameters. The interaction 
term g\ <fi x 2 allows the following tree-level decay process 
to occur in vacuo: — > x + x, if the following mass 
condition is met: m > 2m x . Assuming this condition 
is met, we are led to ask: Will the 0-oscillon decay into 
X? If so, will the growth in x be linear or exponential? 
Otherwise, if m < 2 m x , or if g± = 0, gi 0, we can focus 
on annihilations: (j) + 4> ~ ► X + Xi e ^ c - 

One could approach the issue of multiple fields by 
evolving the full <p, x system under the classical equa- 
tions of motion. In fact by scanning the mass ratio of </) 
and x one can find interesting oscillons involving an in- 
terplay of both fields - one finds that the 2:1 mass ratio 
is of particular importance (as was the case for the SU(2) 
oscillon 4 ). Although this is an interesting topic, here 
we would like to focus on the effects on the 0-oscillon due 
to the introduction of x, initially in its vacuum state. At 
the classical level x will remain zero forever, so this is 
trivial. We will return to the issue of the classical evo- 
lution for non-trivial initial conditions for x in the next 
Section. For now we focus on placing x in its quantum 
vacuum state with a classical background OBC . 

We use the same formalism as we developed in Section 



IV 



We make the replacements 



X, m 



= 1 -> 



in eq. (27), giving the following Heisenberg equation of 



motion for x in fc-space 



Xk 



Wfc Xk 



d d k 
W) 



^(k-k')xfe' =0, 



(67) 



where uj 2 = k 2 + m 2 and $ = gi </> osc (x, t) + g 2 ^„ sc (x, t). 
We write x in terms of its mode functions v q k as before 



Xk(t) = Vfr 



d d q 
(2n) d 



^q Vqk 



(t) + h. 



(68) 



For brevity, lets focus on leading order in e behavior, 
coming from $ g\ 4> e (x) cos(o;i) (although later we 
will also mention the important case of g\ =0 with 
$ = g 2 4> 2 (x)cos 2 (ujt)). This gives the following mode 
function equations 



Vqk+^lVqk+gi COs(wt) 



d d k' 

(2?rp 



s (k-k>,k' =0. (69) 



For small coupling g\ we expect the solutions of the 
mode equations to be small deformation of planes waves. 
To capture this, lets expand the mode functions in powers 
of gi (analogously to the earlier expansion in Section [VJ) 



Vqk 



(0) , (1) , 2 (2) . 

V + 9l Vqk +9lVqk + 



(70) 



At leading order 0(g®), we have v 
.on is th< 

«2(*) = 



(0) 
qk 



,2„,(0) 



v k = 0, whose 



o ~~~~~ - qk 1 ' 

desired solution is the unperturbed mode functions 



(2 7 r) d <5 d (q-k), 



(71) 



as earlier. At next order we have the following forced 
oscillator equation 



7> (1) +^ 2 7) (1) - 
V qk +UJ k V qk - 



e(q-k) 



e -i(ui+u q )t e i(ui-u q )t 



The solution with boundary conditions v i\\o) 



(72) 



••^(0) = is 

"qk 



qk 



(*) = 



E (q-k) 



S(u),U q ,Uk) 



(73) 



where S was defined in cq. (44). Hence we obtain the 



same expressions for the energy density Su and 5E as in 
eqs. (45 46) upon making the replacements 



_Ab_ 
3!2 3 



&(q-k) 



.91 



e(q-k), 



(74) 



Evaluating this we find results that are qualitatively sim- 
ilar to before: the x field is produced in the core of the 
oscillon and spreads out, and the energy grows linearly 
in time. 

Following through a similar calculation to before in the 
small e limit, by identifying |0 e (q — k)| 2 as proportional 
to a (5-function, i.e., 



d d xu e (x) 



(27r) d <5 d (q-k), (75) 



allows us to evaluate the decay rate at leading order in 
<?!, which is connected to the decay process <j) — > x + X- 
We also carry though the calculation at leading order in 
gi (including next order in g\ ) , which is connected to the 
annihilation process (j> + </>—> x + X- We find the results 



fJi 



2 hTT d / 2+1 k d z 2 



2 2 r(f)(27r) 



(76) 



, , . (9i-92) 2 hK d/2+1 kt- 2 U d xu 2 Ax) / . 



r(|)(27r) d / d d xu e {x) 



In (|76|) the kinematical requir eme nt on fc rad is u>k 
k 2 



mr x = (u/2) and in (77) the requirement is 
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J k — n, rad -r in, — lu 2 . So these decays only occur for 
sufficiently light \- Here 1^(1 </> — > 2%) coincides with 
the perturbative decay rate of <j> and 1^(2 — > 2 \) is a 
generalization of the annihilation rate r = n (av) applied 
to a gas of non-relativistic particles with variable density. 
As an application, we expect this result to be relevant to 
the bosonic SU(2) oscillon [4]. Since it exists at the 2:1 
mass ratio (m H = 2m w ), it prevents Higgs decaying into 
W-bosons, but it should allow the quantum mechanical 
annihilation of Higgs into relativistic W-bosons. 



VII. EXPONENTIAL VS LINEAR GROWTH 

So far we have worked to leading order in the couplings 
and e, this has resulted in a constant decay rate of the 
oscillon into quanta of <j) in the single field case or quanta 
of x m t ne two field case. These quanta have an energy 
that grows linearly in time, see eq. (52). However, one 



may question whether this result applies at late times. 
In the case of a homogeneous background pump field, 
it is always the case that the growth is exponential at 
late times if the daughter field is bosonic. This is due to 
a build up in the occupancy number in certain A-modes, 
leading to rapid growth for fields satisfying Bose-Einstein 
statistics. 



A. Floquet Analysis 



particular consider the following pair of (matrix) initial 
conditions 



(<) 
(it) 



Wqk(0) 
W qk (0) 



o, 



Zqk(0) 



0. 

Sqk- 



(81) 
(82) 



Let N be the number of q, k values in our discretization. 
Lets organize this information into a 2N x 2N matrix 
M(t), whose upper left quadrant is w qk with IC (i), upper 
right quadrant is z qk with IC (i), lower left quadrant is 
w q k with IC (ii), and lower right quadrant is z qk with IC 
(ii). So initially we have M(0) = 1-2N,2N- Numerically, we 
evolve this through one period, giving the matrix M(T). 
After n oscillations, we have M(nT) = M(T) n . Hence 
the matrix M(T) controls the behavior of the system. To 
obtain the result for the initial conditions of (30 31 1 we 



multiply M(nT) onto the following diagonal matrix 



vs. 



kq 



o 



o 



-vs~ kq 



(83) 



The existence of exponential growth is governed by the 
eigenvalues of M(T), with some corresponding eigenvec- 
tor {wk, Zk}. Following the standard Floquet theory, we 
write the eigenvalues as exp(/iT), where /j, are the Floquet 
exponents, which are in general complex. Note that al- 
though we have phrased this in the context of solving the 
quantum problem, it is also a classical stability analysis, 
as we mentioned at the end of Section llVl 



Such exponential growth is obtained by fully solving 
the mode functions non-perturbatively [31]. Since the 
background is periodic, the mode functions satisfy a form 
of Hill's equation, albeit with infinitely many coupled os- 
cillators due to the spatial structure of the oscillon. There 
is a large literature on resonance from homogeneous back- 
grounds, particularly relevant to inflation, but rarely are 
inhomogeneous backgrounds studied as we do here. 

The late time behavior is controlled by Floquet ex- 
ponents [i. To make this precise, consider the discrete 
(matrix) version of the mode function equations 

Wqk = z qk (78) 

Zqk = ^2 Q kk ' ( 79 ) 
k' 

where 

Q kk ,(t) = -Jl5 kk , - i$(k- k'.t) (80) 

is a matrix in fc-space, with period T = 2tt/uj. Here we 
have labelled the mode functions w instead of v, since we 
will impose slightly different initial conditions on w. In 



B. Results 

We have carried out the numerical analysis for differ- 
ent models, but would like to report on the results for the 
\g\4>X 2 theory with $ ps g\ (j> £ (x) cos(wt) in d = 1. The 
numerical results for the maximum value of the real part 
of ^i max as a function of g± is given in Fig. [4] (top panel, red 
curve). As the figure reveals, there is a critical value of 
the coupling gl 0.2y/\ governing the existence of expo- 
nential growth. For gi > g\ exponential growth occurs, 
but for gi < gl it does not; in the latter regime all Floquet 
exponents are imaginary. This is reflected in the evolu- 
tion of the energy SE(t), which we have plotted in Fig. E] 
(lower panel). The lower (orange) curve is for g\ = g\[2 
and the upper (green) curve is for g\ = 3(^/2. Hence 
for sufficiently small couplings the perturbative analysis 
is correct - the growth is indeed linear at late times, but 
for moderate to large couplings the perturbative analysis 
breaks down - the growth is exponential at late times. 

How can we understand this behavior? The answer 
resides in examining the structure of the corresponding 
eigenvector. For g\ > gl, we have numerically computed 
the eigenvector {w k , Zfc} max corresponding to the max- 
imum Floquet exponent. It is useful to represent this 
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FIG. 4: Top panel: the maximum value of the real part of the 
Floquet exponent /i max as a function of coupling gi/y \ for 
e = 0.05. Lower (red) curve is for an oscillon pump /i max and 
higher (blue-dashed) curve is for a homogeneous pump ^h,m ax - 
Lower panel: total energy output 8E (xA/<?i) as a function of 
time for <?i/vA = 0.1 (orange) and g\/y\ = 0.3 (green). 

vector in position space, where it is some wave-packet. 
The real part of the (unnormalized) eigenvector is plotted 
in blue in Fig. [5] (top panel). In red we have also indi- 
cated the shape of the oscillon. We find that the shape 
of the wave-packet which carries /i max is approximately 
described by the function 

XmaxW ~ <l>e{%) COS(fc rad X) (84) 

where k lad is the wavenumber we identified in the per- 
turbative analysis (uj\ = k 2 + m 2 , — w/2). Notice that 
the shape of this is independent of the coupling g\. Such 
a result cannot make sense at arbitrarily small values of 



g\. At sufficiently small g\ all eigenvectors of the Flo- 
quet matrix M(T) should be small deformations of plane 
waves, since we are then almost solving the Klein-Gordon 
equation. In particular, this means there should not be 
any localized wave-packet eigenvectors of the Floquet ma- 
trix. If the eigenvectors are spatially delocalized, they 
cannot grow exponentially, since there is nothing avail- 
able to pump the wave at large distances from the oscil- 
lon. This explains why all fi are imaginary at sufficiently 
small <?i. Conversely, at sufficiently large g\, some solu- 
tions can exist that are 0(1) deviations from plane-waves, 
namely the wave-packet of Fig. [5j Clearly then it is inap- 
plicable to treat this as a small perturbation from a plane 
wave. This explains why exponential growth can occur 
at sufficiently large coupling. 

With this understanding, let's postdict the critical 
value of the coupling in this model. If we ignore the 
spatial structure and treat the background oscillon as 
homogeneous with amplitude </> e (0), then %'s mode func- 
tions satisfy a Mathieu equation, whose properties are 
well known (let's call their Floquet exponents /i h ). In 
the regime of narrow resonance, the first instability band 
(connected to <f) ~ ► X + X decay) has a maximum growth 
rate /v max ~ <7i</> £ (0)/2 (plotted as the blue-dashed curve 
in Fig. [I| . On the other hand, the spatial structure of the 
oscillon means that modes that are produced in the core 
of the oscillon will try to "escape" at a rate set by the in- 
verse of the oscillon's width. Let's define an escape rate 
as ficsp — l/(2i? OBC ), where R osc is the oscillon's radius. 
The critical g\ can be estimated by the condition 

<,„ ax * A*-c (85) 

To achieve exponential growth, we require /v max ^ Mesc 
in order for there to be sufficient time for growth to occur 
in the core of the oscillon before escaping, allowing Bose- 
Einstein statistics to be effective. Using 4> e (0) = 4e/vA 
and 1/R OBO ~ e, gives g\ s» vA/4, in good agreement with 
the full numerical result. This reasoning can be extended 
to other scenarios. For g\ = 0, we can focus on anni- 
hilation driven by ^g2(j) 2 X 2 - I n this case, study of the 
Mathieu equation reveals /x h max « g 2 </'e(0) 2 /8, leading to 
i?2 ~ A/(4e). Since A and g\^ are independent parame- 
ters, the regime gi, 2 > g\ 2 is allowed (and easily satisfied 
for XH <C 1, as required for massive oscillons). 

If the parameters are in the regime of exponential 
growth it is interesting to note that substantial paramet- 
ric resonance can occur from an inhomogeneous clump of 
energy established by oscillons. This is a form of paramet- 
ric resonance - explosive energy transfer from a localized 
clump to a daughter field. Of course this cannot continue 
indefinitely, since the oscillon has only a finite amount of 
energy to transfer. Using the initial conditions of Fig. [5] 
(top panel) we have evolved the full coupled {0, x} S Y S ~ 
tcm under the classical equations of motion. We find 



13 




-50 

Position x 



150 



0.15 




FIG. 5: Top panel: The \ wave-packet which exhibits expo- 
nential growth for g\ > gl (blue) and the </>-oscillon (red). 
Lower panel: The fields at t = 80 after classical evolution, 
with gi = 0.8\/A, e = 0.05, and m x = 0.3. 



that exponential growth in \ occurs initially and even- 
tually this results in the destruction of the 0-oscillon as 
seen in Figure [5] (lower panel) . 

Finally, we return to the single field oscillon. In the 
A5 4> 5 model we originally discussed, the resulting (gener- 
alized) Mathieu equation reveals /^ h max ~ A5 </> £ (0) 3 , lead- 



of small amplitude oscillons against exponential growth 
of short wavelength perturbations. This fact was previ- 
ously only seen empirically. It is quite interesting that at 
sufficiently small amplitude, or couplings, the oscillon is 
stable against exponential growth in perturbations and 
yet it still has modes that grow linearly with time. This 
occurs in the limit of degenerate eigenvalues of M(T). 
These modes seem relatively rare and harmless classi- 
cally, but they must be integrated over in the quantum 
theory, resulting in steady decay. 



VIII. COLLAPSE INSTABILITIES 

Some oscillons are unstable to perturbations with 
wavelengths comparable to the size of the oscillon (k — 
(9(e)). Although this is not the focus of our paper, we 
would like to briefly discuss this phenomenon for com- 
pleteness. These instabilities are so prominent that they 
often appear in classical simulations starting from initial 
conditions away from the "perfect oscillon profile" , given 
by the expansion eq. (10). 



ing to A5 ~ A 3/2 /e 2 . However, in our expansion we im- tex and a single contact term from the (j) 4 vertex 



Using the numerical method of the previous Section 
we have found that for k = 0(e) there are exponentially 
growing modes in d = 3, but not in d = 1, for the —A 4> 4 
theory. Since the corresponding wavelengths are compa- 
rable to the size of the oscillon, these instabilities are eas- 
ily seen in numerical simulations starting from arbitrary 
initial conditions; an example is displayed in Figure [6] for 
the Vi ~ —X(f> 4 potential in d = 3. In the simulation, we 
find that the field is growing in the core of the oscillon and 
we also find that it is spatially collapsing. The existence 
of this instability is well known in the literature. 

To feel more confident that this is the correct behavior 
of the quantum theory, let us turn now to address this 
problem in a different approximation. At some level we 
should consider the oscillon as a collection of ^-particles 
in the full quantum theory [34]. Since these "collapse" 
or "self-focussing" type of instabilities occur at small 
wavenumbers k — 0(e), perhaps we can interpret this 
as resulting from the interaction of non-relativistic par- 
ticles. The leading interaction is 2(f) —> 2(f> scattering. 

Lets consider Vj = h A3 3 + |f A4 cf> 4 + At tree- level 

this scattering process occurs due to 4 diagrams: s, t, 
and u-channels generated by 2 insertions of the 3 ver- 

The 



plicitly assumed A5 was 0(1) w.r.t e and therefore we 
should never enter the regime A5 > Ajjj. Similarly, con- 
sider the classic — A cf> 4 model. The resulting (generalized) 
Mathieu equation reveals /i h . max ~ A 2 e (O) 4 ~ e 4 . Com- 
paring this to fi esc ~ e we see that it is impossible to 
obtain exponential growth for small e. (This can change 
for the wide flat-top oscillons of Ref. [H].) As far as 
we are aware, this is the first explanation of the stability 



matrix element is easily computed in the non-relativistic 
limit iM(2 4> — > 2 </>) = — i (|A§ — A 4 ) . This can be recast 
in position space as a two particle potential V(ri, r 2 ) by 
taking the inverse Fourier transform and multiplying by 
h 2 /A due to our normalization convention. This gives 



V(ri,r 2 ) = — 



A 2 - A 4 



5 d {vi-v 2 ). 



(86) 
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Time t 

FIG. 6: Field at center of oscillon over time for Vi = ~4> 4 /4l 
in d = 3. We choose e = 0.05 and set up initial conditions with 
4>(0, p) = 1.03 e 4>i(0, p), i.e., an initial profile 3% higher than 
the preferred oscillon profile; this clearly causes an instability 
by t ~ 150. 



(The s, t, and u-channels actually produce a type of 
Yukawa potential if the full relativistic M. is used, but 
this shouldn't modify our conclusions.) This is a short 
range force and is attractive if and only if §A§ — A4 > 0, 
which is a condition for the existence of oscillons. (This 
condition also emerges in the classical small e expansion) . 
In the non-relativistic regime, the wavefunction for 
particles making up an oscillon will be governed by the 
Schrodinger equation 



N4, 

E v ? + E v ^ r *) I M = En * ^ ( 8? ) 



2m0 



i<3 



with the above potential V. It is known that a 5-function 
potential permits unbounded solutions for d > 2 (d = 2 is 
marginal, being only a logarithmic divergence, but d > 3 
is a power law). A localized gas of ^-particles, such as 
the oscillon, would be expected to be unstable to radial 
collapse under these conditions. Hence, we expect d > 2 
to be unstable (with d = 2 marginal) . For d = 3 a collapse 
time of ~ l/(XHe ) is naively expected. 

There is evidence in the literature [TU [22] that oscillon 
stability is controlled by the derivative of the oscillon's 
mass M aBC w.r.t amplitude <p a . If the derivative is pos- 
itive (negative), then the oscillon appears to be stable 
(unstable) to collapse. Since the canonical small ampli- 
tude oscillon satisfies M osc ~ 4>a~ d > we obtain a consistent 
result. However, this leading order behavior must break 
down at some point for sufficiently large amplitudes and 
can be affected by the inclusion of higher order terms in 



the potential. In fact it is known that beyond a criti- 
cal amplitude no collapse instability exists in 3-d [THl 122"] 
(also see the Q-ball literature [35H37] ). Collapse in 3-d 
is also absent in other field theories, such as the SU(2) 
sector of the standard model [H |S] and any model where 
4> t ~ e 2 instead of the canonical </> e ~ e (e.g., see |18j). 

In summary, the collapse of an oscillon is highly model 
dependent and can be avoided by operating in the appro- 
priate number of dimensions, parameter space, or field 
theory. However, the radiation we computed in the pre- 
vious sections is unavoidable. 



IX. CONCLUSIONS 

We have found that even though an oscillon can have 
a mass that is much greater than the mass of the in- 
dividual quanta, the classical decay can be very differ- 
ent to the quantum decay (this point does not appear to 
have been appreciated in the literature, for instance see 
the concluding sections of Refs. [5, 6 ). The radiation of 
both classical and quantum oscillons can be understood 
in terms of forced oscillator equations, see eqs. (15 41). 



We derived the frequency and wavenumber of the out- 
going radiation, which were both 0(1) in natural units. 
Since a classical oscillon has a spread which is 0(1/ e) in 
position space, it has a spread which is 0(e) in fc-space. 
Its Fourier modes are therefore exponentially small at the 
radiating wavenumber and hence such radiation is expo- 
nentially suppressed. In the quantum theory, there sim- 
ply cannot be modes whose amplitudes are exponentially 
suppressed. Instead, zero-point fluctuations ensure that 
all modes have at least 0{H) amplitude-squared due to 
the uncertainty principle. 

We derived a formula for the quantum lifetime of an 
oscillon ~ l/(Xhe p ). The power p is model dependent: 
p = 4 in the 3 + . . . theory (or — <f> A + 5 + . . .), and 
p = 6 in the — </> 4 + (f> e + . . . theory. Through a Floquet 
analysis, we explained why the growth of perturbations of 
small amplitude oscillons is linear in time, as opposed to 
exponential. The dimcnsionless A h controls the magni- 
tude of the decay rate, as it should for a leading order in 
h analysis. For example, the Standard Model Higgs po- 
tential has XH — (m H /w E w) 2 ~ 0.1 (m H /100 GeV) 2 and 
so this is not very small. On the other hand, the effec- 
tive h of the QCD axion potential is A h ~ (A QCD // a ) 4 ~ 
10 -48 (10 10 GeV/,/ a ) 4 and so oscillons formed from ax- 
ions, called "axitons" in [5], are governed by classical de- 
cay (ignoring coupling to other fields). 

We further considered the fate of an oscillon that is 
coupled to a second scalar x and found it to either decay 
or annihilate with a growth in \ that can be exponen- 
tially fast, depending on parameters. Since oscillons may 
form substantially in the early universe [101 lllj this may 
give rise to interesting phenomenology. At the very least, 
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it presents a plausible cosmological scenario in which a 
parametric pump field exists that is qualitatively differ- 
ent to the homogeneous oscillations of the inflaton dur- 
ing p/reheating. This is a form of parametric resonance: 
explosive transfer of energy from a localized clump into 
bosonic daughter fields. We expect decay into fermions to 
be quite different (for discussion in the context of Q-balls, 
see [35] )■ This may have some cosmological relevance. 

It appears that if a field has a perturbative decay chan- 
nel, then the oscillon will eventually decay through it. 
This is important because we expect most fields in na- 
ture to be perturbatively unstable, including the infla- 
ton, p/reheating fields, Higgs, and most fields beyond the 
standard model. A good exception is dark matter. This 
conclusion may seem surprising given that the oscillon is a 
bound state of particles with a finite binding energy |34j . 
However, oscillons are formed from fields whose particle 
number is not conserved. One could imagine a situation 
in which is only slightly greater than 2m x , and in 
this case the oscillon's binding energy may prevent direct 
decays into %'s, but this requires fine tuning and will not 
forbid 20— !-2xor30— > 2 <fi or 4 <fi — >2<j> annihilations. 

We conclude that in many scenarios an individual os- 
cillon's lifetime will be shorter than the age of the uni- 
verse at the time of production (this may prevent indi- 



vidual oscillons from having cosmological significance in 
such cases). Exceptions include the grand unified theory 
era [11] , inflation [13] , and axitons produced at the QCD 
phase transition [8 . An interesting question for further 
study is whether oscillons can form and then decay, and 
then form again repeatedly, like subcritical bubbles in 
hot water. It is not implausible that such a process could 
continue over long time scales for cosmic temperatures of 
order the field's mass; similar to the production and dis- 
appearance of unstable particles in a relativistic plasma. 
This may modify cosmological thermalization. 
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